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We investigate theoretically the properties of a trapped 
gas of fermionic atoms in both the normal and the superfluid 
phases. Our analysis, which accounts for the shell structure of 
the normal phase spectrum, identifies two observables which 
are sensitive to the presence of the superfluid: the response 
of the gas to a modulation of the trapping frequency, and 
the heat capacity. Our results are discussed in the context of 
experiments on trapped Fermi gases. 

The observation of Bose-Einstein condensation in sev- 
eral atomic systems [|l| has recently sparked increasing 
interest in trapped fermionic atoms. These systems of- 
fer the prospect of a Bardeen-Cooper-SchriefFer (BCS) 
transition to a superfluid phase at low temperatures 
T < Tc- By trapping the atoms in two hyperfine states, 
the phase transition temperature Tc should be experi- 
mentally accessible [H,|| , and several experimental groups 
are presently working to achieve this transition [0,p|. 
However, as only a few percent of the atoms are likely to 
participate in Cooper pairing 0, it is not obvious how 
the transition could be observed in these dilute systems. 
Recently, it has been proposed that the propagation and 
scattering of light should be significantly altered by the 
presense of Cooper pairs 0,§. Since the quasiparticles 
(QP) with energies near the Fermi chemical potential, 
/ip, are those most affected by the Cooper pairing, can- 
didate observables for the detection of the BCS transition 
should be sought from phenomena sensitive to this low- 
energy region of the QP spectrum. 

In this paper, we consider two such observables: the 
response of the gas to a "shaking" of the trap, as first 
suggested by Baranov |^ ; and the heat capacity. For low 
T, both of these observables are dominated by contri- 
butions from the low-energy spectrum. By presenting a 
complete calculation of the properties of the trapped gas 
in both the normal and superfluid phases, which accounts 
exactly for the quantization of the single-particle energy 
levels, we are able to predict if these two observables are 
suitable to detect the presence of Cooper pairing. Our 
analysis should have direct relevance to the ongoing ex- 
periments on trapped Fermi gases. 

We consider a gas of fermionic atoms of mass m, con- 
fined by a potential C/o(r), with an equal number of 
atoms TVo- in each of two hyperfine states, jtr = ±). Two 
fermions in the same internal state a must have odd rela- 
tive orbital angular momentum (minimally p-wave) , and 
at low temperatures the centrifugal barrier suppresses 
their mutual interaction |1C| ]. Thus, we suppose the in- 
teraction to be effective only between atoms in different 
hyperfine states and to be dominated by the s-wave con- 
tribution. As the interplay between the discrete nature of 



the normal phase spectrum and the Cooper pairing is cru- 
cial for the quantities considered in this paper, we need 
a theory which can describe the interaction and Cooper 
pairing of atoms residing in different discrete trap levels. 
This precludes the use of a simple Thomas-Fermi treat- 
ment pl| . A theory appropriate for the present paper has 
recently been presented ||l^. It uses a zero-range pseu- 
dopotential 11^ to model the interaction between atoms 
in two different hyperfine states; this is appropriate when 
the scattering length for binary atomic collisions, a, has 
a larger magnitude than the effective range of the inter- 
action, Te, and when kp\a\ <C 1, where k^ = y^m/xp/^ i^ 
the Fermi wavevector. The generalized mean field theory 
derived from this approach yields the eigenvalue prob- 
lem f2i: 



E,,u^{v) = [7^0 + W{v)]u^{y) + A(r)«^(r) 
E^v,,{v) = -[no + W{r)]vr,{r) + A(r)u^(r). 



(1) 



Here Ho = ^|^V^ + Uo{r) — /ip is the single-particle 
Hamiltonian; W(r) = g('i/;J(r)^CT(r)) is the Hartree po- 
tential, where ^'^(r) is the atom field operator for com- 
ponent a at position r, which obeys the usual fermion 
anticommutation relations. The coupling constant is 
g = Anafi jvfi and the pairing field, A(R), is defined 
by 

A(R) ^ -g lim 9,[r(^>+(R + \)^-{^ - \))\. (2) 

Our definition of the pairing field differs from that of- 
ten employed in weak-coupling BCS theory W%- The 
main advantage of the definition given above is that it 
eliminates the ultraviolet divergence present in the usual 
weak-coupling theory |12(| . The elementary quasiparti- 
cles (QPs) with excitation energies £^^ are described by 
the Bogoliubov wave functions Ur^ij) and i'r;(r)- 

We solve the Bogoliubov-de Gennes (BdG) equa- 
tions (0) for the case of an isotropic harmonic po- 
tential, C/o(^) = TOa;^r^/2, using a self-consistent nu- 
merical procedure outlined elsewhere [|2|. In the ab- 
sence of the pairing field, the QPs exhibit a discrete 
spectrum of energies E',,, with the index 77 designat- 
ing a triple of quantmn numbers {n,l,m), where l,m 
are the usual angular momentum quantum numbers and 
n is an index of radial excitation. In the presence 
of the pairing field, the self-consistent solution to the 
BdG equations with the lowest free energy is spheri- 
cally symmetric. Thus, the Bogoliubov wavefunctions 
are given by u,,(r) = r~^Uni{r)Yim{0,4>) and Vn(r) = 
T^^Vni{r)YimiO,4>)^ where the YJm are the usual spheri- 
cal harmonics, and with n, I, m being implicitly indexed 



by r]. The paring field, which is a scalar operator un- 
der rotations, couples a normal phase QP with {n,l,m) 
to one with (n' , I, ~m). Due to the spherical symmetry, 
in taking sums over states needed to obtain the results 
of this paper, we can replace sums over m by factors of 
(2/ + 1). 

We now calculate the response of the gas to a harmonic 
time-dependent perturbation of the trapping potential, 
AH{t), of the form 

An{t) ^ Asin((::>i) ^ fd^r^mLj^r^^jl{r)i;„{r), (3) 

where A is a small parameter. We expand the field op- 
erators in terms of the Bogoliubov wave functions and 
the QP operators in the usual way p4| , and by applying 
Fermi's golden rule to obtain the linear response R{uj) of 
the gas to the perturbation, Eq.(||), we obtain: 



/•oc 

R{w) oc2 Y^ (2^ + 1)1 / dr{i 



[UnlUn'l - VnlVn'ljr 
n>n',l 

X ifn'l - fnl)5{hLb + En'l - Enl) + 

V" (2/ + 1)1 / dr{u,iiv,i'i + VniUn'iy^f 
^^ 1 Jo 

X (1 - fnl - fn'l)Sihi:j - Enl - E„,i), 



2|2 



(4) 



where /„; = (exp /?£;„; + 1) \ /3 = l/fceT, and fee is 
Boltzmann's factor. The physical interpretation of the 
two terms in Eq.fl) is straightforward: The first term 
describes the excitation of a QP due to the perturbation, 
whereas the second term describes the creation of two 
QPs. This latter process does not violate particle con- 
servation, since the QPs in general are mixtures of real 
particles and holes. The response of the gas should be 
observable as density fluctuations of the trapped gas. As 
we have assumed a spherical symmetric perturbation, the 
transitions all have Al = 0. A generalization to pertur- 
bations with arbitrary angular momentum / is straight- 
forward. In the non-interacting limit, Eq.(0) reduces to 
a sum of delta functions S{CJ — 2nuj) with n = 0,1,2 .. .. 
We now solve the BdG equations self-consistently and 
then calculate the response of the gas to a "shaking" 
of the trap from Eq.(|4|). In Fig.|l|, we show a typ- 
ical plot of the response R{oJ) for various values of 
T' = k^T/hia. In this example, we have chosen the 
parameters g/{hujlf^) = —0.8 and /zp = 51.5huj, where 
Ih = {h/mioY^^ is the characteristic length of the ground- 
state harmonic oscillator wavefunction. With the value 
of a = — 2160ao, appropriate to ®Li |1^, these parameters 
correspond to Nc^ ~ 3.8 x 10'' atoms of each spin state 
in a trap with frequency v — uj/2'k ~ 520Hz; a value of 
Tc ~ 5.6huj/kB = 140nK for the transition temperature 
is obtained by linearizing Eq.(|l|) ||l^. Fig.|l| shows the 
response for T' = 0, 3.95 and 4.55, where the gas is in 
the superfluid phase, and for T' = 6 > T^, where the gas 
is in the normal phase. For comparison, we also plot the 



T = response, assuming the gas is in the normal phase. 
Each delta function in Eq. (|j) representing a i ^ oo reso- 
nance, is smoothed out to a frequency range of ^ tj/10 to 
model the finite frequency resolution of the appropriate 
experiment. 

We now discuss these results, considering first the re- 
sponse for the normal phase. The resonance peaks for 
T = and T' = 6 are relatively narrow on the scale of 
Lo. This is perhaps surprising, as one might expect the 
Hartree field to wash out the shell-structure of the QP 
spectrum in the normal phase Pq |. To understand this, 
we plot in Fig.g the lowest QP energies Ey^ for the gas 
in the normal phase at T = 0. To simplify the plot, 
we include only even values of I. The QP energies with 
odd I behave in a completely analogous way. All ener- 
gies are positive; negative normal-phase particle energies 
are simply holes {u^^{r) = 0) with positive energy in this 
representation. For T = 0, only the 5(?icj — E^j — Ejj') 
term in Eq.(|4|) is non-zero. The thick vertical arrow in 
Fig.S indicates a typical transition: creation of a hole 
with energy Eh, and a particle with energy Ep, yielding 
huj = Eh -\- Ep ^ 2.2huj. The analysis of this normal- 
phase spectrum is basically the same as the one presented 
in Ref . |l6| . A key result of the present paper is the find- 
ing that, although the Hartree field has introduced a sig- 
nificant dispersion of the QP energies as a function of I, 
the dispersion is almost the same for each band. Hence, 
for AZ = 0, the difference of energies between two par- 
ticle bands (or the sum of energies of a particle and a 
hole band, as in Fig.g) varies much less with I than the 
energies themselves, which results in a relatively narrow 
resonance peak. The resonance for T = is sharper than 
for T' = 6. This is because for T — Q, only the energy 
bands immediately around /ip contribute to the response 
due to the Fermi exclusion principle, whereas for higher 
T, there are transitions between several bands that yield 
slightly different transition energies. 

We now consider the response when the gas is in the 
superfluid phase. By comparing the result for T' = & 
and T' = 4.55 in Fig.|], we see that when the gas en- 
ters the superfluid phase, there is a significant broaden- 
ing of the resonance line. This is due to the fact that 
Cooper pairing starts to mix particles with holes, and 
the QP spectrum is altered. This is depicted in FigJ^, 
which shows the lowest even-Z QP levels for T' = 4.55 for 
both superfluid and normal phases. When the energies 
of particles and holes are almost degenerate in the nor- 
mal phase {I ~ 26 in Fig.||), the pairing strongly mixes 
these two states. This leads to the usual avoided cross- 
ing and the QP-spectrum is changed significantly. The 
strong mixing yields the broadening of the resonance line 
depicted in Fig|l|. There are now transitions with signif- 
icantly lower energies than in the normal phase. Such 
a transition, which contributes to the 5{'hLJ + E^y — E^,) 
term in Eq.(0) with Ej^ — Er,' ~ 0.8?iw, is indicated by 
the vertical arrow in Fig.|^. We also note that the effect 
of the pairing decreases with increasing I. This is sim- 
ply because the centrifugal potential "pushes" the high I 



states into the region where the order parameter becomes 
very smaU. The inset in Fig.|l| shows A(r) and |M^(r)|. 
For T' = 4.55, the pairing only takes place around the 
center of the cloud, and QP states which have a small 
amplitude in this region are unaffected. 

For T -^ Tc, all the low lying QP states are strongly 
influenced by the pairing. From the inset in FigJ^, we 
see that Cooper pairing now takes place over the entire 
trapped cloud. The low energy QP spectrum for T = 
plotted in Fig.g is qualitatively different from the normal 
phase spectrum. The low energy QP wave functions are 
centered between the regions where the pairing field and 
the trapping potential are significant. These "in-gap" 
states, which were first discussed by Baranov j|], depend 
strongly upon the strength of pairing. As T decreases and 
A(r) increases, their energy increases. The response of 
the gas is completely dominated by these states for T <C 
Tc- The broad peak for T = in FigjI] comes from the 
S{huj — En — En') term in Eq.(^) with E^ — E^' being the 
lowest energy for a given I. It reflects excitations of the 
kind 7^„7^_^|$o) where |<i>o) is the ground state and 7^^ 
creates a QP with quantum numbers 77 in hyperfine state 
a. Hence, for T ^ Tc the response of the gas is a broad 
peak coming from excitations of the lowest QP band. The 
resonance peak should be centered around an increasing 
frequency as T is lowered, since A(r) increases. This is 
confirmed in Fig.nl where a broad peak has emerged in 
the response for T' = 3.95, the peak being centered at a 
lower frequency than for T = 0. 

The qualitative behavior of the response of the gas 
described above depends on the fact that the resonance 
peaks are relatively well-defined in the normal phase. We 
have performed a number of calculations varying both 
the coupling strength and the number of atoms trapped. 
For experimentally realistic parameters, it turns out that 
the Hartree field does not wash out the resonance peaks 
in the normal phase. We therefore believe the analysis 
above should be valid for typical experimental conditions. 

The low-T heat capacity is another observable which 
probes the low lying QP spectrum. The usual way to 
measure the energy of a trapped gas is to turn off the 
trapping potential and then deduce the velocity distri- 
bution from the expanding cloud ||I7| . As the trapping 
potential is turned off non-adiabatically, the energy ob- 
served is really Etot — Epot where Etot is the total energy 
of the trapped gas and Epot — X^cr /^'^''Po-('")"^'^^^^/2 
with Pair) = {tpl{r)ipa-{r))- Thus, the most appropriate 
definition of the heat capacity for the present purpose is 
Cjv = dxiEtot — Epot)\N„- Solving the BdG equations 
with varying T, we can calculate Etot{T) and Epot{T) 
and therefore Cn ■ The total energy of the gas given the 



solution of Eq. (|lD is 



I 



+ Y.(^ - f,) J d'rv^{r){n, - EnWir). (5) 



As an example, we plot in Fig.^ cn = Cn/'^Nct both 
for the normal and the superfluid phases. The number of 
particles in the trap is held constant at 2Na- — 24860 and 
g / (fiLolf^ — —1, corresponding to a trapping frequency of 
820Hz for ^Li. The critical temperture for this set of pa- 
rameters is UbTc ~ 4.5?laj. From Fig.^, we see that the 
heat capacity is suppressed in the superffuid phase for low 
T . This is because the pairing removes the gapless ex- 
citations present in the normal phase. These excitations 
have equal particle and hole character, and are there- 
fore strongly infiuenced by the pairing as noted earlier. 
Therefore, the heat capacity is exponentially suppressed 
by a factor ^ exp(/3A), where A is the gap in the QP 
spectrum coming from the Cooper pairing. However, the 
suppression is only significant for T <^Tc, where all angu- 
lar momentum states are affected by the pairing and the 
QP spectrum is truly gapped (compare Fig.y and FigJ^). 
Since the system is finite and the superfiuid correlations 
occur gradually starting in the center of the trap at Tc^ 
and then continuously extending outwards as T decreases 
(see inset of Fig.p, there is no discontinuity in the heat 
capacity at T = Tc, in contrast to the case of an infinite 
homogeneous system Q. It should be noted that it is 
important that the normal phase spectrum is approxi- 
mately gapless on the scale of A such that the normal 
phase Cn behaves linearly with T for T < Tc as depicted 
in Fig.H. Otherwise, the heat capacity is exponentially 
suppressed even in the normal phase, and one will not 
observe a significant change when the gas becomes su- 
perffuid. Fortunately, for realistic values of g and N„^ it 
turns out that the Hartree field indeed makes the normal 
phase QP spectrum essentially gapless |16| . 

In conclusion, we have presented a detailed analysis of 
two possible ways of detecting the predicted BCS phase 
transition for a trapped gas of fermionic atoms. The on- 
set of Cooper pairing infiuences significantly the response 
of the gas to modulation of the trapping frequency. For 
T > Tc, the response has a relatively sharp peak, and 
the width of the peak should narrow as T is lowered. 
Then as T = Tc is reached, one should observe a sig- 
nificant broadening of the response peak as the Cooper 
pairing starts to affect the low lying QP spectrum. For 
T <^Tc, the low lying QP states are qualitatively differ- 
ent from the normal phase states, and the response of the 
gas to the shaking is predicted to be a broad peak coming 
from the lowest QP band. The center of the peak should 
move to increasing frequencies as the pairing increases 
for decreasing T. Also, one should be able to detect the 
phase transition by looking at the low T heat capacity. It 
should be exponentially suppressed for T <SiTc refiecting 
the gapped nature of the QP spectrum due to Cooper 
pairing. However, a measurement of the heat capacity is 
destructive as one has to release the trap and it requires 
several repetitions of the trapping experiment. Also, the 
suppression of Cn^ is only significant for T <^ Tc- We 
therefore estimate that it is a less direct way of detecting 
the transition than by looking at the response to a modu- 
lation of the trapping frequency. The analysis presented 



here should be quaHtatively correct for a non-spherical 
symmetric trap as well, although the actual calculations 
would be more cumbersome in this case. 



Rb; M. H. Anderson et al, Science 269, 198 (1995); Na: 
K. B. Davis et al, Phys. Rev. Lett. 75, 3969 (1995); 
Li: C. C. Bradley et al, ibid. 78, 985 (1997); H: D. 
G. Fried et al, Phys. Rev. Lett. 81, 3811 (1999). See 
http://amo.phy.gasou.edu for references to other ex- 



[9 

[lo; 
[11 

[12; 

[13 

[14 

[is; 
[16; 

[17 



porimcnts. 

H. T. C. Stoof, M. Houbiers, C. A. Saclcett, and R. G. 
Hulet, Phys. Rev. Lett. 76, 10 (1996) 
M. Houbiers, R. Ferwerda, H. T. C. Stoof, W. L 
McAlexander, C. A. Sackett, and R. G. Hulet, Phys. Rev. 
A 56, 4864 (1997) 

D. S. Jin, Centennial Meeting of the American Physical 
Society, March 20-26, 1999, Atlanta, Georgia, Vol. 44, 
No. 1, FB16 8 

M.-O. Mewes, G. Ferrari, F. Schreck, C. Salomon, Cen- 
tennial Meeting of the American Physical Society, March 
20-26, 1999, Atlanta, Georgia, Vol. 44, No. 1, FB16 7 
M. Tin- 

kham. Introduction to Superconductivity (McGraw-Hill, 
New York, 1975) 



J. Ruostekoski, cond-mat/9902324 

Z. Wang, C . A. Sackett, and R. G. Hulet, cond- 



mat/9902350 



M. A. Baranov, |cond-mat 9801142| 

B. DeMarco and D. S. Jin, |cond-mat/9807406 

D. A. Butts and D. S. Rokhsar, Phys. Rev. A 55, 4346 
(1997) 

G. M. Bruun, Y. Castin, R. Dum, and K. Burnett, Euro. 

J. Phys. (in the press) 

K. Huang, Statistical Mechanics (Wiley, New York, 

1987). 

P. G. de Gennes, Superconductivity of Metals and Alloys 

(Addison- Wesley, New York, 1966). 

E. R. I. Abraham, W. I. McAlexander, J. M. Gerton, R. 
G. Hulet, Phys. Rev. A 55, R3299 (1997) 

G. M. Bruun, and K. Burnett, Phys. Rev. A 58, 2427 

(1998) 

J. R. Ensher, D. S. Jin, M. R. Matthews, C. E. Wieman, 

and E. A. CorneU, Phys. Rev. Lett. 77, 4984 (1996) 



DC 



T=0, Normal Phase 




-- 


~X 


\ / 


5 20 




\ 
\ 




"g 




\ 


T'=6 ll 


Sio 


-^ 


x^--^.. 


\ll 





5 10 I J 1 T=0 




'"" r 1 1\ \ 




T =3.95 \l\ \ \ \ 




'' y^ ■'■ ^'** ^^ 


T'=4.55 


>:;..:(1. J%^, \ 



00/ CO 

FIG. 1. The response K of the gas as a function of the 
modulation frequency iI) for various temperatures. The in- 
set shows the Hartree field |W^(r)| (dot-dashed) and A(r) for 
T' = (solid) and T' = 4.55 (dashed) in units of ht^. 
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FIG. 2. The lowest QP energies in units of ?itj at T = for 

the normal phase (x) and the superfluid phase (Solid line). 
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FIG. 3. The lowest QP energies in units of fiuj at T' = 4.55 
for the normal phase (x) and the superfluid phase (Solid line). 




FIG. 4. The heat capacity in units of fcs for the normal 
(dashed) and the superfluid (sohd) phase. 



